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Abstract 

Surface and curvature properties of asymmetric semi-infinite nuclear matter are 
studied to beyond the proton drip. Using the semiclassical extended Thomas-Fermi 
method, the calculations are performed in the non-relativistic and relativistic mean 
field theories (Skyrme forces and non- linear a—uj parametrizations). First we discuss 
the bulk equilibrium between the nuclear and drip phases. Next we analyze the asym- 
metric surface as a function of the bulk neutron excess. We examine local quantities 
related to the density profiles and, for two definitions of the bulk reference energy, the 
surface and curvature energy coefficients. The calculation of the curvature energy is 
carefully treated. The sensitivity of the nuclear surface to the relativistic effects is in- 
vestigated. Mass formulae useful for arbitrary neutron excess are discussed, and their 
limit at small asymmetries is compared with the liquid droplet model mass formula. 



PACS: 21.60.-n, 21.10.Dr, 21.65.+f 
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1 Introduction 



In several problems of nuclear physics and astrophysics the surface and curvature prop- 
erties of nuclei play a crucial role. This is the case, for instance, of barrier heights and 
saddle-point configurations in nuclear fission, or of fragment distributions in heavy-ion colli- 
sions. In astrophysical applications they are important for equilibrium sizes, electron capture 
rates and level densities used in describing neutron stars and supernovae. However, the pro- 
ton fraction of the nuclei involved in these problems is rather different. For terrestrial nuclei 
the proton fraction is around 0.4-0.5, whereas the astrophysical problems demand a smaller 
proton concentration and, furthermore, to consider nuclei surrounded by an external gas of 
drip neutrons. Therefore, there is a strong motivation for studying the surface and curva- 
ture properties of nuclei in the whole range of asymmetries, between nuclear and neutron 
matter. 

Within a context related to the liquid droplet model (LDM) and the leptodermous ex- 
pansion of a finite nucleus fl]], the surface and curvature energy coefficients can be extracted 
from semi-infinite nuclear matter. This schematic geometry avoids undesired shell. Coulomb 
and finite-size effects, which makes it very appropriate to study surface properties. Many 
of the calculations have been carried out using the semiclassical extended Thomas-Fermi 
(ETF) method 0,0,^ together with Skyrme forces 0. For the case of symmetric semi- 
infinite matter, the reader can consult e.g. Refs. Although the surface energy can 
be calculated in a fully quantal way, only one part of the curvature energy can be obtained 
from a quantum-mechanical calculation @,|^. This is an important reason to compute the 
surface and curvature coefficients by means of a semiclassical approach. 

As far as asymmetric semi-infinite nuclear matter is concerned, we may point out two 
publications of particular relevance for the present work. The surface energy of multi- 
component systems up to the proton drip was investigated by Myers et al. |]TU[ in the 
Thomas-Fermi (TF) approximation with the Seyler-Blanchard interaction. Almost at the 
same time, Kolehmainen et al. |Tl[] performed an analysis of the surface and curvature effects 
of neutron-rich nuclei, including drip neutrons, with the ETF method and Skyrme forces. 

In recent years the relativistic treatment of the nuclear many-body problem has been a 
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subject of growing interest |T2|,|T3|,|T^,|T5[ . The relativistic theory, already in the mean field 
(Hartree) approximation, automatically includes the spin-orbit force and the finite range and 
density dependence of the nuclear interaction. All these effects are very important for the 
surface properties of a nuclear system [^. The phenomenological a—uj model of Walecka 
and its extensions |lT2| , p!7| have become very popular in relativistic nuclear calculations. They 
have been widely employed in the Hartree approach for describing ground-state properties 
of nuclei [|rB|,|r^,^,^ as well as symmetric semi-infinite nuclear matter [ P^j2^J23[ . 

Until recently the semiclassical approach to the relativistic mean field theory had been 
settled at the TF level only [Il2| , p!7| , p4| , p5| , p6| . However, in the last years some amount 
of work has been addressed to develop the relativistic extended Thomas-Fermi (RETF) 
method p7|^^p0| , that includes gradient corrections of order in the particle density 
and the scalar field. This new semiclassical approach has already been applied to the study 



of finite nuclei and of the nuclear surface in the symmetric case P^ JS 

In comparison with quantal calculations, the inclusion of the gradient corrections to the 



TF method improves the surface properties (energy and profile) ||31|j3^j35| . Moreover, the 
quality of the RETF results is less dependent on the parameters of the relativistic interaction 
than in the TF approximation. 

The asymmetric matter has also been investigated within the relativistic theory. For 
example, the infinite system has been analyzed in the Dirac-Hartree-Fock approach [Q, and 
surface properties at low proton concentrations have been computed with the TF method 
37 1 . Miiller and Dreizler |^| have calculated relativistic hot asymmetric nuclear matter 
and nuclei, and in Ref. Miiller and Serot have carefully dealt with the subject of phase 
transitions in asymmetric nuclear matter at finite temperature. Other works on relativistic 
asymmetric systems have been concerned, e.g., with neutron star properties or collective 
modes and response functions 

The main aim of this paper is to analyze, from a semiclassical point of view, the surface 
and curvature properties of asymmetric semi-infinite nuclear matter in the non-relativistic 
and relativistic frameworks. For this purpose we will use the ETF and RETF methods 
respectively. Our study covers the range of proton concentrations from symmetric nuclear 
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matter to beyond the proton drip. We will put special emphasis on the calculation of the 
curvature energy, that can be fully obtained only within the semiclassical approach and in 
which some confusion exists [p|,p|,|35|. 

The paper is organized as follows. Section 2 is devoted to the study of uniform asym- 
metric nuclear matter with arbitrary proton concentration, including drip neutrons and 
protons. In the third section we derive the equations for the interface surface and curvature 
properties in the general case with drip particles. Numerical investigations are presented 
and discussed in Section 4. The summary and conclusions are given in the last section. 
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2 The bulk equilibrium 



Before proceeding to study the surface properties of nuclei in coexistence with drip 
nucleons, we start by the simpler model of two infinite pieces of asymmetric nuclear matter 
in equilibrium. The uniform phase with the higher density and a given proton concentration 
represents a nucleus that is in equilibrium with a surrounding gas of drip nucleons. The 
latter is represented by the lower density phase, with another proton concentration. This 
approach neglects the interface effects and reflects the situation of the equilibrium densities 
very far from the interface region in either direction. 

The asymmetric nuclear matter problem including drip particles has a notable simili- 
tude with the problem of the coexistence between nuclei and evaporated nucleons at finite 



temperature [pq,H3,H^,44]. Actually, the bulk equilibrium of asymmetric matter at finite 



temperature was investigated time ago within the astrophysical context |45,46|, and more 



recently in Refs. |^9[^. In this section we perform a similar analysis at T = based on 
the relativistic non-linear a—uo model, with care over the effects of the proton concentra- 
tion. However, our discussion is rather general and can be applied to the non-relativistic 
treatment of cold asymmetric nuclear matter as well. 

First we consider an infinite system made up of two components (neutrons and protons) 
in a single phase. For describing asymmetric matter it is convenient to switch variables from 
the neutron and proton densities pn and Pp to the total density p and the relative neutron 
excess 5: 

P = pn + Pp (2.1) 

^^Pn^_ (2.2) 
P 

If we call S{p,6) the energy density of the system, then the neutron and proton chemical 
potentials and the pressure are 

_ dS{p,6) _ d^dp_ d^d6__dS_ dS 1 - 6 
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Manipulating these expressions, one gets that the pressure is related with the neutron and 
proton chemical potentials through 



P = PnPn + PpPp - £{P, . 



(2.6) 



At saturation P = and the preceding equation is, actually, a generalization of the 
Hugenholtz-Van Hove theorem ||48| , |49[| . The incompressibility at saturation of asymmet- 
ric nuclear matter with relative neutron excess 6 reads 



K{p,,t,S) = 9p^-^ 



rfp2 1^ p ^ 



(2.7) 



Psat 



where psat is the saturation density that conforms Eq. (|2.6| ) with P = 0. The expression 
of the energy density S for the relativistic model and for Skyrme forces can be found in 
Appendix A. There, £ is given including the semiclassical gradient corrections of order fi^ 
which, of course, vanish in the uniform infinite geometry. 

The asymmetric nuclear matter in a single phase is not stable at all densities and S 
values. At zero temperature, the necessary and sufficient conditions for the stability of a 
binary system can be expressed through the inequalities p9| , ^ , |50 
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In addition, the pressure P must be positive. A positive incompressibility K guarantees 
mechanical stability, and the condition on the derivatives of the chemical potential ensures 
diffusive stability. Violation of the stability criteria signals phase separation. The system 
ceases to exist as a phase alone and splits into separate phases in equilibrium. However, 
instability is only a sufficient condition and phase separation may take place when the 
single-phase system is stable (in general, metastable). 

Now we address the problem of the coexistence of two phases in a two-component system. 
One of the phases represents an asymmetric uniform nuclear system with density po and 
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relative neutron excess 60 . The other phase corresponds to a uniform phase of drip particles 
characterized by pd and 6^- Equation ( p.6|) is still valid for each separate phase, but at 
equilibrium the pressure is no longer zero. In the general case with drip neutrons and 
protons, the two phases can coexist when the following conditions are fulfilled: 

A^nO = /ind (2.10) 

yWpO = f^pd (2.11) 

^^0 = ^^.1, (2.12) 

which in view of Eq. (|2.6|) result in 

/in(PnO - Pnd) + l^piPpO " Ppd) = Sq - Sd ■ (2-13) 



Equations (|2.1(]| )- (|2.12|) allow one to get po, pd and 6d for a fixed value of 60 in the nuclear 



phase. Notice that the solution of the above equations does not correspond to the saturation 
conditions at 60 and 6d- As a consequence, the energy of the nuclear and drip phases in 
coexistence is not a minimum. 

The numerical results of this section will be discussed for the NLl parametrization of 
the relativistic non-linear a—uj model ||20|. Other cr—uj parameter sets or non-relativistic 



Skyrme forces qualitatively show the same trends. The saturation properties of NLl are 
given in Table 1, together with those of other interactions we will utilize later in the calcu- 
lation of surface properties. 

To study the different regimes of phase coexistence. Fig. 1 displays the neutron and 
proton chemical potentials against the relative neutron excess for NLl. At small values of 
So there are no drip particles and the nuclear matter fulfills the usual saturation condition. 
In this first stage the chemical potentials are negative. While for neutrons the chemical 
potential is an increasing function of 60, it decreases steadily for protons. The point where 
the neutron chemical potential vanishes marks the onset of the neutron drip. At this moment 
the coexistence between the nuclear and drip phases, the latter containing only neutrons 
(^d = 1), starts to be possible. In this second regime, that for NLl covers the range 
0.24 < Sq < 0.64 approximately, the neutron chemical potential (/ino = Pnd) is positive. 



Since there are no protons in the drip matter yet, the proton chemical potentials fipo and 
/ipd of the nuclear and drip phases are different. Both fi^Q and yUpd are negative, being 
lowered by raising the neutron excess. Due to the fact that fipd is larger and decreases at 
a faster rate than yUpo, they match at some large Sq, which for NLl happens at 6q ^ 0.64. 
This is just the point where protons start to drip. The last two rows in Table 1 display the 
values of 60 at the neutron and proton drip points (labelled and 5pd). 

For larger 60 the drip phase contains a small proton concentration {6^ < !)• The two 
phases disappear when the solution of Eqs. (|2.10|) -( p.l2|) yields the same density and neutron 
excess in the nuclear and drip regions. This solution is just the critical point characterized 
by pc and 6c, see Eq. ( |2.14] ) below. It corresponds to the maximum of fino{6) and to the 
minimum of fipo{6). For 6 larger than 6c the dashed lines in Fig. 1 represent the neutron 
and proton chemical potentials of the drip matter only. Tracing horizontal lines in the fi — 6 
plane one can read, at the left and right hand sides of the critical value 6c, the relative 
neutron excess of the nuclear and drip phases in coexistence beyond the proton drip point. 

Figure 2 shows the nuclear and drip densities, po and pd, as a function of the square 
of the neutron excess for NLl. For small values of 6q the nuclear density po has a linear 
behaviour with 6q, as expected from the droplet model The density of the nuclear (drip) 
phase falls (climbs) with increasing 6^. Both densities coincide at the critical 6c- The dashed 
line of Fig. 2 represents the density of the drip region for 5^ > 5^. As in Fig. 1, by drawing 
horizontal lines it serves to read the equilibrium values of 6^ and 6^ when protons have 
migrated to the drip phase. 

The equation of state of NLl at several values of 6 is displayed by solid lines in Fig. 3. The 
dot-dashed line defines the boundary of the coexistence region. The long-dashed and short- 
dashed lines are the diffusive and mechanical stability lines respectively. We have joined 
by dashed horizontal lines a few examples of points on the coexistence curve, belonging to 
the nuclear and drip phases, that can be in equilibrium. The common maximum of the 
coexistence and diffusive stability lines corresponds to the critical point {6c, Pc) determined 

by 
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8 



The line of mechanical stability lies in the negative pressure region. Therefore, it does not 
play any role in the separation of the system in two phases, that is completely ruled by the 
diffusive instability. The zone between the coexistence and diffusive stability lines is the 
metastable region where the binary system may remain in a single phase or undergo phase 
separation. 

Consider a uniform nuclear system with density p and neutron excess 5 under a pressure 
P. If the representative point in the P — p diagram of Fig. 3 lies outside the coexistence 
region, the system is in a single phase. If the system is expanded along a line of constant 5 
it will remain in a single phase up to crossing the coexistence line, where phase separation 
can occur. Denoting by ^ the mass fraction of the nuclear phase (5o, po) in coexistence with 
the incipient drip phase {5^, p^), with mass fraction 1 — ^, the following relationships are to 
be satisfied at that point: 

5 = ^So + il-O^d (2.15) 
1 £ 1 - £ 

- = ^ + — ^. (2.16) 

P Po Pd 

In Figs. 4 and 5 we plot in the Pn — P and Pp — p planes the lines of coexistence (solid), 
diffusive stabihty (long dash) and mechanical stabihty (short dash) for NLl. The chemical 
potentials and densities at the right and left hand sides of the critical point (the maximum 
and the minimum of the coexistence line in the Pn — p and pp — p planes) correspond to 
the nuclear and drip phases respectively. The dashed straight lines joining points of the 
coexistence curve indicate some possible values of the densities and chemical potentials in 
the phase equilibrium. Notice that in the Pp — p plane some of these straight lines are not 
horizontal. In such cases the drip phase does not contain protons {Sd = 1) and the proton 
chemical potentials of both phases are different, see Fig. 1. When the proton drip point is 
reached these lines turn horizontal, indicating that Ppo = Pp^- 
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3 Surface properties in two-component systems 

On the basis of asymmetric semi-infinite nuclear matter, in this section we shall study 
the surface effects of a neutron-rich nucleus immersed in a gas of drip particles. This regime 
corresponds to the physical situation found in neutron stars, at densities between 4 x 10^^ 
g/cm^ and normal nuclear values (2.7 x 10^^ g/cm^) First of all it is necessary to 

determine the density profile from which the surface and curvature energies can be obtained. 
To do that one considers a semi-infinite slab with a plane interface separating a mixture 
of protons and neutrons to the left and a gas of drip particles to the right. The axis 
perpendicular to the interface is taken to be the z axis. The density profile is sketched in 
Fig. 6. Observe that the relative neutron excess b depends on the z coordinate. When z 
goes to minus infinity, the neutron and proton densities approach the values that in uniform 
nuclear matter are in equilibrium with a uniform gas of drip particles. Therefore, for a given 
bulk neutron excess 5o = ^(~oo), the values of p(— cxd), p(oo) and b{po) are just the po? Pa 
and (5d solutions to Eqs. (|210|)-(|2A^). 

The problem of finding the density profile at a certain value of 5o when drip particles 
appear is, in a sense, formally equivalent to solving the problem of a warm nucleus. In 



that case, as it was pointed out by Bonche et al. the Hartree-Fock equations at some 
temperature and chemical potential possess two different solutions. One of them represents 
the nuclear system in equilibrium with its evaporated nucleons. The other one belongs to 
the evaporated nucleon gas alone. It is then natural to define extensive magnitudes char- 
acterizing the nucleus as the difference between their value in the nuclear-plus-gas solution 
and their value in the gas solution. The same idea has been applied in the semiclassical 



scheme for the Euler-Lagrange equations at finite temperature . 

To compute semiclassically the proton and neutron densities in asymmetric semi-infinite 
nuclear matter one has to minimize the total energy per unit area. This has to be done 
with the constraint of conservation of the number of neutrons and protons with respect to 
arbitrary variations of the densities. According to the above discussion, the energy, neutron 
and proton densities entering the constrained energy of the (semi-infinite) nucleus will be 
taken as the difference between the ones of the whole system (nucleus plus drip particles). 
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and those of the uniform gas of drip particles (indicated by a d subscript). For instance we 
will write the local energy density representing the nucleus, i.e. the subtracted system, in 
the form £{z) — As z ^ — oo we have £{z) — » £q, while as z — » cxo it is £{z) £^. 
Due to the absence of Coulomb forces in the semi-infinite problem, the drip particles are 
described by plane waves and behave as a dilute uniform phase. 
Then, we write the constrained energy per unit area as 

= dz{£{z) -£d- fin[pn{z) - pnd]- fip[pp{z) - ppd]} ■ (3.1) 

For the whole system, the proton and neutron densities are obtained from the following 
coupled Euler-Lagrange equations: 

^-/in = 0, — l^-/ip = 0. (3.2) 

In the relativistic framework one has three additional equations, originating from the vari- 
ations with respect to the fields 0, V and R associated with the a, uj and p mesons: 

IT-' IT-- 

We recall the reader that the ETF and RETF expressions for £ are detailed in Appendix 
A. The drip particles obey analogous equations: 

/in = , /ip = , (3.4) 

C>Pnd ^Ppd 

plus in the relativistic case 

o0d oVd oRa 



The local equations ( |3.2| ) and ( ^.31 ) are self-consistently solved by numerical iteration. 
As explained in Ref. [28|, we employ the imaginary time-step method to get the densities 
from Eq. ( |3.2| ) and Gaussian elimination to get the meson fields from Eq. ( |3.3| ). Given 
the asymptotic boundary conditions to be imposed stem from the coexistence equations 



(|2.1CI| )- (|2.12|) . Equations ( p.4|) and ( p.5|) are satisfied by the low-density phase solution of 



( |2.1CI| )-( ^.12| ), which provides the limiting behaviour as ^ oo. On the other hand, the 
high-density solution of (p.lCI|) - (p.l2|) represents the system at z = — oo. 
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Once the density profiles are known, the next step is to get the surface and curvature 
energy coefficients. The surface tension (i.e. the surface energy per unit area of a flat surface) 
may be written as [^,1C,11] 



/oo 
dz[£{z) - Sd- Srcliz)], (3.6) 
-oo 

where Sj-ef{z) is a reference energy density whose integral represents in some way the bulk 
contribution. Myers et al. [|10| noted that there are two possibilities to define £re{{z)- 

The first definition corresponds to a reference energy that represents the energy a nucleus 
would have if its nucleons would belong to infinite nuclear matter. This reference energy 
density is to be written as 

= §^ W) - Pd] , (3.7) 

PO — Pd 



which we will call the e-definition following Ref. |10]. In this case the surface tension reads 



(Jo 



dz \£{z) -g^- ^^-^ [p{z) -pA. (3.8) 

-oo [ Po — Pd J 



The second definition of the reference energy will be called the Gibbs definition. Instead 
of the bulk energy per particle, it introduces the neutron and proton chemical potentials 
and the pressure associated with the bulk. The Gibbs reference energy for the whole system 
is 

E^^, = /i„iV + /ipZ - PV, (3.9) 

while for the drip phase we have 

= PndiVd + Ppd^d - P^V. (3.10) 

The meaning of E'^ becomes evident if one writes the thermodynamic relation 

dE = p^dN + ppdZ - PdV , (3.11) 

where dE is the energy necessary to remove dN neutrons, dZ protons and reduce the volume 
by an amount dV. Thus, E^^ may be interpreted as a reference disassembly energy [|T0| . 



12 



With the help of the equihbrium conditions ( p.l0| )-( |2.12| ) between the nuclear and drip 



phases, the reference energy for the subtracted system reads 

Ki-KfA = /in(iv-iVd) + /ip(z-Zd) 

/oo 
dz {Hn [Pn{z) - Pnd] + /^p [pp{z) - Ppd] } • (3.12) 
-oo 

From this equation we can identify a Gibbs reference energy density: 

Cef(^) = Pn [pn{z) - pnd] + Pp [pp{z) - ppd] ■ (3.13) 

Substituting Eq. ( p.l3| ) into Eq. (|3.6| ), the Gibbs surface tension becomes 

/oo 
dz{£{z) -Ea- pApA^) - Pnd] - Pp\pp{A - Ppd]} • (3.14) 
-oo 

As a side remark we note that £^^^^£(—00) = E'^^,^{—oo) = So — S^ and that £^^^^£(00) = £^^'^£(00) = 
0. 

Comparing with Eq. ( p.l| ) one concludes that the surface tension that must be minimized 
for obtaining the density profile is the Gibbs surface tension (|3.14| ). Attempts to minimize 



the surface tension cXe, Eq. ( |3.8| ), result in the non-conservation of neutrons and protons 
Nevertheless, cXe can be calculated from the variational densities arising from the 



minimization of cr^. Of course, in the case of symmetric semi-infinite nuclear matter, ae and 
cr^ coincide owing to the Hugenholtz-Van Hove theorem. The surface energy coefficients 
associated with and cr^ are given by 

Es,e = 47rroVe (3.15) 
E,,^ = 47rroV^, (3.16) 

where tq is the nuclear radius constant: 

^nrl{po-pd) = l. (3.17) 

Another useful quantity for our study is the so-called surface location zq. It is defined 

as 
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_ JZ.zp'{z)dz 
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where the prime denotes a derivative with respect to z. Similar expressions hold for the 
neutron and proton surface locations (2;on and Zop)- It is easy to show that 



zoAp = 2;onApn + ^OpApp 



(3.19) 



with Ap = po — pd, Apn = Pno — Pnd and App = ppo — Ppd- After some algebra one obtains 
an expression relating the two definitions of the surface tension: 

, Apn App 



(A'n — A'p)(^On — ^Op)" 



Ap 



(3.20) 



In semi-infinite nuclear matter the curvature energy coefficient is given by p|, ^ , ^ 

/oo 
dz{z-z^)\E{z)-EA-E,,,{z)\ 
-oo 

d 



/oo (-1 
^dz — [S{z)-S^~S,,i{z)] 



(3.21) 



ft=0 



where k is the curvature {k, = 2/R for a sphere of radius R). The two contributions to the 
curvature energy coefficient in Eq. ( |3.21|) are called geometrical (-Ef''°) and dynamical (E^^^) 
respectively. The geometrical contribution only involves the variation of the surface energy 
density S{z) — — Sj-cfiz) across the surface parallel to the z axis. The dynamical part 
comes from the change of the surface energy density by curvature when the plane surface is 
infinitesimally bent. 

In order to establish a connection between E^^° and E^^^ it is helpful to define the surface 
width of the density: 

2 _ IZojz - zoyp'{z)dz 
j_^p'{z)dz 

plus similar quantities and bp for the neutron and proton densities. One finds that b, b^ 
and 6p satisfy 



6" (Ap)' = bl Ap^Ap + bl AppAp + {zon - z^^f Ap^App . 



Making some calculations it can be shown that 



2Ap 



^n~\^ (^On — Zqp) 



2 Apn - App 

Ap 



(3.23) 



(3.24) 
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The dynamical part of the curvature energy requires paying a special attention. In 
general, the surface energy density £{z) — — Sj-cf{z) depends on the curvature k in two 
different manners. First, S{z) may have an explicit dependence on the laplacian operator 
A which in the limit R ^ oo reads as (P /dz^ + nd/dz. Second, the particle densities (and 
the meson fields in the relativistic case) carry an implicit dependence on the curvature k. 
By construction the energy density entering Eq. ( p.21| ) is free from any explicit dependence 
on the Laplace operator (see Appendix A). Actually, this dependence has been removed 



by partial integrations in the semiclassical functionals P8| , |35| . Consequently, in the present 
calculation, the sole contribution we may have to the dynamical curvature energy comes from 
the implicit curvature dependence of the nuclear densities and fields. For non-relativistic 
interactions we find 



Et^^ = SnroYl H dz \ — [S{z) - - S.^z)] ^ 
„rr„j-oo I 0,, dn 

(3.25) 



K = 



In the relativistic case one has a similar expression but including the variations with respect 
to the meson fields |]9|,|35[]. 

If we insert in Eq. ( |3.25| ) the Gibbs reference energy E'^^^i^z)^ Eq. ( |3.13D , the prefactors of 
the derivatives with respect to /t vanish by virtue of the variational equations obeyed by the 
self-consistent densities (and meson fields), cf. (p.2|)-( p75D . Therefore, in our calculation the 
Gibbs dynamical curvature energy is zero [E'^y^ = 0) and the full Gibbs curvature energy 
coincides with its geometrical part: Ec^f_i = Ef^". 

This is not the situation when one takes the e-definition of the reference energy £'^^^{z), 
Eq. ( p.7| ). In this case one finds a non-zero dynamical curvature energy -E^e"' because the 
prefactors of the derivatives on k in Eq. (|3.25|) do not vanish, (-^c e" — only for symmetric 
matter with 5o = 0.) Notice that in the relativistic framework there are no contributions 
from the meson fields to E^^f (as long as we employ the functionals of Appendix A), since 
they do not enter in the reference energy. 

E^^^ can be recast in the form 

roo Q 



/oo cj 
Jz — [£{z)-£,-£U^)] 



K=0 
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+ 87rro 



(3.26) 



K=0 



The first term is just E^^^, tliat vanishes according to the foregoing discussion. As a 
consequence, the dynamical part of the e-curvature energy can be expressed as 



<?=n,p 



^0 — 



oo Q 
-oo OK 



(3.27) 



K = 



Po - Pd. 

In general the evaluation of dpg/dn in Eq. ( |3.27] ) is not a trivial matter, e.g. it cannot 
be calculated in a fully quantal way . Fortunately, this problem can be solved within the 
semiclassical formalism 1^,0. The reason is that the semiclassical expression of the particle 
density derived from the Wigner-Kirkwood h expansion of the density matrix contains 
Laplace operators, and thus a dependence on the curvature k. From the Wigner-Kirkwood 
density, following Refs. and the dynamical curvature energy E^^J^ can be calculated 
in the ETF and RETF approaches. This is summarized in Appendix B. 
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4 Discussion of Results 



To illustrate our results, we have chosen the Skyrme forces SkM* and Sill as 
representative interactions for the non-relativistic case, and the non-linear parameter sets 
NLl ||2^, NL2 and NL-SH ||5^ in the relativistic model. Table 1 collects the saturation 
properties of nuclear matter for these interactions, plus the droplet model parameters J and 
L: 

2p d6^ 



5=0, p=p, 

3p d fS{p,6)\ 



(4.1) 



<5=0,p=pnn 



2 dpdS"^ \ p 

with pnm the saturation density of symmetric nuclear matter. Also given in Table 1 are the 
values and 5pd of the bulk neutron excess at the neutron and proton drip points. 

Although SkM* and Sill resemble in the energy per particle, effective mass m^/m 
(~ 0.8) and bulk symmetry energy J (~ 30 MeV), these forces differ mainly in the incom- 
pressibility modulus K (217 MeV for SkM* and 355 MeV for Sill). In the relativistic case, 
NLl, NL2 and NL-SH also disagree in the value of K: NLl is similar to SkM*, NL-SH to 
SIII, and NL2 has the largest K (~ 400 MeV). The NLl and NL2 sets have a higher bulk 
symmetry energy J (~ 45 MeV) than the Skyrme forces, whereas in NL-SH the value of 
J is relatively close to SkM*. Concerning the effective mass, in the relativistic sets it is 
small as compared with the Skyrme forces. For NLl and NL-SH m^/m is similar (~ 0.6), 
while for NL2 it is slightly larger. In any case, notice that the effective mass has a different 



origin in the relativistic than in the non-relativistic model |£8|. All the forces are able to 
give a reasonably good description of finite nuclei in spite of their differences. In particular, 
the relativistic set NL-SH is very well suited for nuclei near the neutron drip line that 
cannot be described so well with other a — uj sets or Skyrme forces. 



4.1 The surface properties 

Figure 7 displays the neutron and proton local density profiles obtained from the solution 



to the Euler-Lagrange equations ( 3.2|) -( p75D for NLl and SkM*. They are plotted for several 
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values of the relative neutron excess in nuclear matter 60 . The separation between vertical 
bars is the surface thickness t, defined as the 90% to 10% fall-off distance. At the top of the 
figure we show the result for 60 = 0, the symmetric semi-infinite geometry. The neutron and 
proton profiles coincide and the surface region is roughly centered around z = 0. The second 
part of the figure corresponds to 60 = 0.2 for which the density profiles begin to differ. The 
surface thickness for neutrons and protons is very similar to the 6q = case. However, while 
t is still centered around z = for the neutron profile, for the proton density t is shifted 
to negative values of z owing to the symmetry terms of the interaction. The third part of 
the figure corresponds to 6q = 0.4, when drip neutrons have appeared both in NLl and in 
SkM*. The surface thickness of the neutron and proton distributions has clearly increased. 
The density of the drip phase is considerably higher for NLl than for SkM*, as one would 
expect from the fact that the neutron drip starts earlier for NLl (Table 1). Finally, at the 
bottom of Fig. 7 the profiles are plotted for a high Sq beyond the proton drip point. The 
surface thickness becomes quite large because the interior density and the density of the 
drip phase are very close. 

Figures 8 and 9 offer a more detailed analysis of the surface. They display the surface 
thickness t and the surface width b, Eq. (p.23|) , for neutrons and protons calculated with 
the relativistic sets NLl and NL2 and with the Skyrme forces SkM* and SIII. One can 
see that the slope of the t and b curves as a function of 60 is steeper for neutrons than for 
protons. After the neutron drip point, one observes the appearance of a relative maximum 
(whose height strongly depends on the force) in the surface thickness and width of the 
neutron density distribution. The presence of this maximum has already been reported in 
earlier literature 0]. From a qualitative point of view it can be understood as follows. If 
the bulk symmetry coefficient J is large (as in the a— 00 sets NLl and NL2), it costs the 
system a great amount of energy to produce an asymmetry in the bulk, thus favouring the 
ejection of neutrons to the surface region. However, for a given 60 the neutron density in the 
drip region is fixed by the coexistence equations (|2.10| )- (|2.12|) . Therefore, if the neutrons 
pushed out by the symmetry term of the force cannot accommodate in the drip region, they 
concentrate at the surface and contribute to the development of a maximum in t and b. 
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This fact is illustrated in Fig. 10, where the number of neutrons per unit area in the surface 
region N/S is drawn as function of Sq for NL2. We have calculated N/S as the integral 
of the neutron density Pn{z) — within the 90%-10% fall-off distance. One recognizes 
that it is immediately after the neutron drip point {6^^ = 0.222 for NL2) when the greatest 
accumulation of neutrons in the surface takes place. 

Another quantity of interest for inspecting the surface of the asymmetric system is the 
neutron skin thickness: 

9 = Zon - , (4.2) 
where zq^ and zop are the neutron and proton surface locations defined through Eq. ( p.l8|) . 



Figure 11 depicts the change of against 6q for our representative forces. For small values 
of 6q the neutron skin thickness behaves linearly, as predicted by the liquid droplet model 
(LDM) i: 

e^^^*. (4.3) 

The value of Q, which measures the stiffness of the system against pulling zq^ and zq^ apart, 
can be obtained from the slope of Q at 6o = 0. At small 6o, forces having a large J/Q ratio 
have also a large neutron skin thickness. As seen from Fig. 11, for larger 6o the neutron 
skin thickness starts to depart from the linear behaviour and develops a maximum in the 
region beyond the neutron drip point. This maximum, which also appears in former analyses 
with non-relativistic forces |rD,rT|, is more pronounced for the relativistic interactions we 
study here. It is related with the accumulation of neutrons at the surface we have discussed 
before. With further increase in 6q the neutron skin decreases, since the inside and the 
outside matter become more and more alike and the surface is washed out. At the critical 
6c one expects O to eventually vanish. 

Figure 12 displays and cr^ against 6q. The surface tension grows with 6q and reaches 
a maximum at some point beyond the neutron drip. Then it falls off to zero as the limiting 
situation of a uniform system is approached. The Gibbs surface tension shows a different 
behaviour and it decreases monotonically from its value in the symmetric case to zero. In 
Figs. 13 and 14 we present our results for the 7e and 7^ curvature energies per unit length 
(7 = Ec/Sirro). These quantities exhibit a similar behaviour as a function of 6q to the surface 
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tensions. NLl presents a visible peak in 7e. The dynamical contribution 7^^" = E'^^^/'^tttq is 
shown by dashed lines in Figs. 13 and 14. Remember that in our calculation the dynamical 
curvature energy only differs from zero in the e-definition for 5q 7^ 0. At high values of 5q 
most of the e-curvature energy comes from the dynamical part. 
For small values of 5q the nuclear droplet model predicts 

i?s,o is the surface energy coefficient for the symmetric case, K is the bulk incompressibility 



and L reads for the LDM coefficient defined in Eq. ( [4.1|) that gives the density dependence 
of the symmetry energy. As we realize from Eqs. ( [4.4| ) and ( [4.5| ) , the symmetry contribution 
to the surface energy in the LDM is positive in the case of the e-definition and negative with 
the Gibbs prescription. This contribution consists of two terms. The main one (9J^/4Q), 
represents the variation of the bulk symmetry energy J when the nucleus increases its 
neutron skin against the resistance provided by the surface stiffness Q. The corrective term 
2Es qL/K describes the change of the volume energy produced by a change in the bulk 
density M. While both terms participate additively in the symmetry contribution to i?s,e, 



they have opposite signs in £'s,^. Since ag and must behave according to Eqs. (^4.41) 
and ([4.5| ) for small Sq, and they must vanish for the uniform system, one can qualitatively 
understand the global trends of Fig. 12. Apart from using Eq. ( [4.3| ) to obtain the surface- 
stiffness coefficient Q, it can be extracted from the slope of the difference -Es,c — E^^^ that 
for small values of 60 behaves as 

E.,e - E,,, = . (4.6) 

Table 2 collects E^, E^ and t computed in the symmetric case {6q = 0) and the coefficient 
Q, for SkM*, SIII, NLl, NL2 and NL-SH. For the sake of comparison. Table 2 shows in 



addition results from Refs. [[1^ and Our values of the surface-stiffness coefficient 

calculated from Eqs. (|4.3| ) and (|4.6| ) are in good agreement between them, though small 
differences arise. In general, the value of Q from Eq. ( |4.3| ) is ~ 0.5 — 1 MeV greater than 
the value from Eq. ( [1.61) . 
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The strength of the peaks that appear in t^, &n and B in Figs. 8, 9 and 11 is, actually, 
related with the value of the surface-stiffness coefficient of the force. Comparing with Table 
2, one can see that forces with small values of Q have the peaks more developed. As 
discussed above, such peaks are connected with the neutrons pushed to the surface. Since 
Q measures the resistance against separating neutrons from protons, it is qualitatively clear 
that forces having a small Q will tend to concentrate more neutrons at the surface. A similar 
dependence on Q is observed in the height of the maxima presented by and 7e in Figs. 
12-14. 

From Table 2 we can see that non-relativistic and relativistic parametrizations are able 
to give comparable surface properties. Nevertheless, for the symmetric case, in Ref. a 
tendency to thinner density profiles was noticed in the relativistic model, with smaller values 
of the thickness t when the relativistic interaction was adjusted to give the surface energy 
of Skyrme forces. Figure 8 suggests that this trend may be reversed at large asymmetries, 
where t (and b as well. Fig. 9) grows faster with 6q for NLl and NL2 than for SkM* and 
SIII. To get more insight, we have performed calculations with the relativistic set named 
RSkl* in Ref. The parameters of RSkl* were fitted to the nuclear matter properties 
of SkM* (including the coefficient J), and the scalar mass was chosen to obtain the surface 
energy of SkM* at 5o = in the ETF approach. For RSkl* we find L = 81.8 MeV and 
Q ^ 25 MeV, whereas for SkM* it is L = 45.8 MeV and Q ^ 39 MeV. Observe from this 
and Tables 1 and 2 an overall tendency of the relativistic sets to have larger L and smaller 
Q than the Skyrme forces. 

Figure 15 summarizes our results for the surface thickness and surface tension with SkM* 
and RSkl*. The relativistic thickness t is smaller for low neutron excess. But it increases 
more rapidly with 6o and becomes clearly greater than for SkM* at high asymmetries, 
meaning that the interface disappears sooner in the relativistic case. The height of the 
relative maximum developed by the surface thickness of the neutron density is larger for 
RSkl* (g ^ 25 MeV) than for SkM* {Q ^ 39 MeV). In Fig. 15 the curves of the Gibbs 
surface tension cr^ for RSkl* and SkM* closely follow each other (cr^ is the minimized 
quantity in the semi- infinite calculation). At low 6o, the surface tension cxe is larger for the 
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set RSkl* than for SkM*, while the contrary happens with cr^. This is so because of the 
values of Q for both interactions, as the LDM equations ( [4 .41) and (|4.5| ) show. Since by 
construction both forces have the same J, the quantity 9J^/AQ is larger for RSkl* and thus 
the initial slope of the relativistic and cr^ is steeper than for SkM*. Altogether, the above 
hints at a difference in the asymmetry dependence of the nuclear surface in the relativistic 
model as compared to conventional Skyrme forces. 

To conclude this discussion we would like to clarify why the relativistic parametrizations 
give larger values of L. This fact can be traced back to the behaviour of the symmetry 
energy with density. From Eq. ( ^4.1| ) and Appendix A, for Skyrme forces it is 

h^kl 1/1 \ 1 /37r2\^/\ ^ ^ ^ 

Jip) = -^-i^o(2+''V''+24 [^2(4 + 5x2)-3tiXi]//3 

with kp = (37rV/2)^/^. The first term comes from the kinetic energy and the other ones from 
the interaction. For the considered forces the term with top gives a positive contribution, 
whereas the terms with p^/^ and p^^° produce a negative contribution. The net result is 
that the function J(p) initially grows with density, but it is progressively bent as the density 
increases. In the relativistic model one has (Appendix A, h = c = 1) 

J(P) = ^ + 1-2P- (4.8) 

The first quantity is the relativistic kinetic contribution, while the isovector p meson is the 
mechanism of the relativistic interaction to introduce additional symmetry energy. In Eq. 
( |4.8| ) there are no negative terms that oppose the growing behaviour of J{p). Since the 
coefficient L is essentially the derivative of J(p) evaluated at pnm, cf. Eq. this explains 

a higher L in the relativistic model. 

4.2 Comparison with previous calculations 

For the Skyrme forces, our results can be contrasted with those of Kolehmainen et al. 
[pT| (KPLT hereafter). Our calculation differs from KPLT in two aspects. On the one 
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hand, our approach is fully variational whereas in KPLT the densities were restricted to 
have the form of a trial Fermi function to a power (plus a constant in the case of drip). 
On the other hand, we have kept the coefficient of the Weizsacker term in the second-order 
energy functional equal to its standard ETF value, namely (3 = 1/36. However, in several 
applications of KPLT it was set to = 1/18 which somehow simulates effects of order h'^ 
[^0. From Table 2 one can see that, for f3 = 1/36, the results obtained in KPLT by means 
of a parametrized density are very close to those we obtain with fully variational densities 
(see also Ref. |^). It is also clear that the use of /3 = 1/18 or /9 = 1/36 has a drastic effect 
on the surface properties. 

For the SkM* and SIII interactions we may compare our results for the surface tensions 
(Te and 0"^, shown in Fig. 12 as a function of 6q, with the corresponding values presented in 
Fig. 7 of KPLT (called cts and Us respectively). The KPLT calculation was performed for 
f3 = 1/18, but we find a fair agreement between both calculations if our results are scaled 
by a factor ctkplt,o/c"o (the ratio at 6^ = of the KPLT value with f3 = 1/18 to ours with 
P = 1/36, that can be read from Table 2). 

From the comparison of Fig. 13 of this work with Fig. 8 of KPLT, we see that the 
agreement is also good for the Gibbs curvature energy per unit length 7^ (called Uc in 
KPLT), when our results are scaled by 7kplt,o/7o (it can be read from Table 2). However, a 
remarkable discrepancy appears when the values for 7e (called cxc in KPLT) are compared. 
After scaling our values with 7kplt,o/7o5 they are larger than those reported in KPLT. The 
disagreement persists if we perform our self-consistent calculation with (3 = 1/18. One must 
note that the authors of KPLT did not include the dynamical curvature energy which comes 
from the implicit curvature dependence of the nuclear densities (see Section 3 and Appendix 
A of KPLT). This is a non- vanishing contribution in the case of 7e at 6q 7^ 0, Eq. ( p.27|) . 

In Table 2 we also compare our surface properties for the relativistic model with those 
reported by Von-Eiff et al. |3^ (VPSW hereafter). The VPSW results were computed at the 
TF level, without the gradient corrections we have in the RETF method. For the surface 
energy coefficient, in the symmetric case, one finds lower values in the RETF than in the 
TF calculation. As discussed in Refs. |p3tp8i^^ , this is a typical feature of the relativistic 
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sets that have a bulk effective mass m*^/m, roughly, below 0.70. From Table 2 we see that 
the surface thickness t turns out to be smaller in RETF than in TP, in accordance with 
the trend of the surface energies. Finally, the surface stiffness Q is larger in the RETF 
calculation than in TF, which expresses a greater rigidity of the nuclear system in RETF 
against the separation of the neutron and proton surfaces. 



4.3 Nuclear matter and the surface in the relativistic model 

Next, we analyze the dependence of the surface on the various magnitudes that charac- 
terize the infinite nuclear matter in the relativistic case. The saturation properties of nuclear 
matter are governed by the meson coupling-to-mass ratios g'^/m'^, g1/m1 and and by 

the non-linear couplings g2 and (73 [p!2| , p]7| (see Appendix A). Contrarily, the surface prop- 



erties extracted from the semi-infinite system depend on the meson coupling constants and 
masses separately. 

We start by obtaining the surface properties with a non-linear a—uo set whose parameters 
we have fitted to these nuclear matter properties: volume energy = — 16 MeV, saturation 
density pnm = 0.16 fm~^, incompressibility K = 200 MeV, effective mass m*^/m = 0.70, 
and bulk symmetry energy J = 30 MeV. The scalar meson mass is rris = 500 MeV. We 
call this set of parameters NLM. Later we recalculate the surface quantities changing one 
of the properties that define the set NLM at each time, with the other properties fixed to 
their initial values. This way we can study the individual effect of each bulk property and 
of the scalar mass rus on the surface. The results are collected in Table 3 for the symmetric 
case 60 = 0, and in Table 4 for 60 = 0.212 (below neutron drip, corresponding to an ideal 
system of ^°*Pb) and 60 = 0.6 (above neutron drip). 

Table 3 also shows the value of the neutron excess in uniform matter when drip neutrons 
appear (5nd) for NLM and the related sets. It can be seen that important changes in 6nd 
are connected with varying the volume and symmetry energies. This is due to the fact that 
the neutron drip is mainly determined by the ratio a^/ J. On the other hand, by decreasing 
the effective mass or the incompressibility, the neutron drip point moves to higher values of 
60. Of course, 6-ad is not changed by variations of mg that are compensated by changes in 
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the scalar coupling constant g^. The proton drip point is less sensitive to the nuclear matter 
properties and appears at 6q ~ 0.7 for all the examined cases, and thus we do not display it 
in Table 3. 

Since the scalar mass sets the range of the attractive scalar interaction, there is a direct 
correlation between rus and the surface properties. A larger rris determines a shorter range 
of the attractive potential. As Tables 3 and 4 show, this leads to a steeper surface and to 
a visible reduction of the surface and curvature energies and of the surface thickness t and 
the neutron skin thickness 0. The surface also is strongly correlated with the value of the 
effective mass of the interaction. Indeed, m^/m plays a prominent role in the majority of 
nuclear structure properties in the relativistic theory, since it is intimately related to the 



vector and scalar fields ||14| , |21| . For the symmetric case and small asymmetries, on increasing 
mj^/m the surface and curvature energy coefficients and t and B show a downward trend, 
as it happened with the scalar mass. The same situation is found at 6q = 0.6 except for 
-Es,^, which turns out to be larger for m*^/m = 0.70 than for m*^/m = 0.55. The effective 
mass provides a measure of the non-local effects, which contribute to make the surface more 
diffuse [0 . A higher effective mass is associated with less non-locality and thus it tends to 
favour a sharper surface. 

A smaller value of the incompressibility K softens the nuclear surface and the thickness 
t augments, while the surface and curvature energies decrease. The change of the neutron 
skin thickness B induced by K is not monotonous with 6o- From Tables 3 and 4 it can be 
checked that the main changes in the ratio Eg/t come from the value of K, as discussed in 
the literature [ p!^j2^ . Increasing K brings about larger Es/t values. This ratio experiences 
a clear reduction for high 6q. For instance, for the set NLM one has Eg/t = 8.8 MeV/fm at 
So = 0, 9.1 (6.0) MeV/fm at 5o = 0.212, and 3.7 (0.4) MeV/fm at 5o = 0.6 (the first number 
corresponds to Eg^e/t and the number in parentheses to Es^^jt). 

Modifying the saturation density pnm or the volume energy has little consequences on 
the surface properties. The analysis of Tables 3 and 4 shows that has a moderate effect 
on the surface and curvature energies, while the surface thickness t is not sensitive to pnm- 
Conversely, the incidence of is more visible on the surface thickness than on the energies. 
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The neutron skin thickness Q is almost unaffected by pnm and a^- 

The coefficient J does not change the surface properties of the symmetric system. The 
surface energy Es^c and the neutron skin B increase with J at 6q = 0.212, but the opposite 
tendency is found at 6q = 0.6. At both values of 60, E^^e and the thickness t become larger 
when J is raised, whereas Eg^^ and i?c,/i decrease. Reading from Table 3 the value of Q 
for NLM and for the set with J = 40 MeV, the behaviour of G, Eg^^, and E'g^^ with J at 



^0 = 0.212 is consistent with the LDM equations 

We can see from Table 3 that the surface- stiffness coefficient Q is raised by larger values 
of rris, m*^/m and J. This behaviour of Q with mg and m*^/m is not surprising because 
the surface becomes stiffer when these quantities increase, and it is harder to separate the 
neutron and proton surfaces. Furthermore, for small 5q at least, it reflects the inverse 
proportionality between Q and the neutron skin B shown by Eq. ( |4.3|) , since G decreases 
with mg and m*^/m (at fixed 5q). Reducing and or a greater K also makes Q grow, 
but the changes are less noticeable. 

In contrast with our finding that Q increases with J when the other properties of the 
force are not changed, in previous literature increasing values of J have been related with 
decreasing values of Q |TT1 , |37|| . Certainly, this is true if the total (bulk-plus-surface) sym- 
metry energy -Esym of the interaction is set to the empirical value fitted by mass formulae, 
that contain a term of the type a^yraiN — ZY /A. Introducing I = {N — Z)/A, in the LDM 
one has |jl| 

E.y^ = a,,^ I^A =[j-[%- ^) A^'^'] I' A , (4-9) 

which is valid for large A and small /. Equation (|4.9| ) shows that if the symmetry energy 
-E'sym of a given interaction is fixed empirically, a large J coefficient must be associated with 
a small Q coefficient and viceversa. On the contrary, if the parameters of the interaction 
are not fitted to reproduce -Egym such constraint between J and Q needs not be satisfied, as 
happens e.g. when the parameters are chosen to give the nuclear matter properties we wish. 

To clarify this point we plot in Figs. 16 and 17, respectively, the surface-stiffness coeffi- 
cient Q and the mass-formula symmetry coefficient agym (calculated for ^''^Pb) as a function 
of J. They are drawn for several a—ui and Skyrme parametrizations; we have added the 
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relativistic set SRK3M5 (J = 23.5 MeV) E^ and the Skyrme force SI' (J = 29.35 MeV) 



ri| , |60| to the interactions aheady considered. It can be seen that forces fulfilling that Q 
decreases as J increases (namely SIII, SI', SkM*, NL-SH and NLl) roughly lie on a curve in 
the Q — J plane (Fig. 16). These same forces lie in a rather narrow region of values of agym 
(22-24 MeV, the empirical region) in the asym — J plane (Fig. 17). There are two forces, 
namely NL2 and SRK3M5, that clearly deviate from the general tendency in the Q — J 
plane. The reason for this anomaly is that these two forces are clearly outside the empirical 
region where the remaining forces lie in the agym — J plane. 

For further insight we have drawn in Figs. 16 and 17 the results of two more forces, 
NL1J4 and NL1J5. They are identical to the set NLl excepting the value of the bulk 
symmetry energy: J = 40 MeV for NL1J4 and J = 50 MeV for NL1J5. The coefficient 
asym for NL1J4 and NL1J5 moves away from the empirical band, but both Q and agym grow 
with J because the remaining nuclear matter properties have been kept constant. Thus, to 
get the right symmetry coefficient agym one must change various properties of the effective 
interaction simultaneously, and then the tendency of smaller Q with higher J is fulfilled. 

4.4 From semi-infinite matter to finite nuclei 

To conclude this section, it may be interesting to compute the energy of large asymmetric 
and uncharged nuclei using the mass-formula coefficients calculated in the semi-infinite 
medium. We can then check the results against those we obtain from a self-consistent. 



semiclassical calculation of spherical finite nuclei p8|j3T| , |59| without Coulomb interaction. 
Before presenting the numbers we pass to discuss some aspects referred to the mass formulae 
we will utilize. 

In writing the energy of a nucleus one has two choices, according to the two definitions 
of the reference energy, Eqs. ( |3.7| ) and ( |3.12| ). One possibility is 

E = Ev,e(5o) A + E,,e(5o) A^'^ + Ec,e(5o) A^'^ , (4.10) 



with an e- volume energy E^ ^ = (^^o ~ (Po ^ Pd), see Eq. (|3.7| ). The other alternative to 



27 



write the energy is 

E = E,,^{So) A + E,,^{So) ^2/3 + E,,^{So)A^/^ , (4.11) 

in terms of a Gibbs volume energy E^^^ = fin{N — N^)/A + fip{Z — Z^)/A, see Eq. ( p.l2| ). 
The surface energy coefficients i?s,e and i?s,^t and the curvature energy coefficients E^^e and 
i?c,/i have been defined in Eqs. (|3.15|) , (|3.16|) and ( p. 211) . 

In order to estabhsh a hnk with the usual mass formulae we shall restrict ourselves to 
values of 6o below the neutron drip point. In this case i?v,e(5o) is just the energy per particle 
in bulk matter S/p, calculated at the saturation density psat for ^o- Likewise, in E^^^{6o) the 
neutron and proton chemical potentials will be evaluated at psat and 6q, with = Z^^ = 0. 

When one deals with finite nuclei it must be taken into account that the overall neutron 
excess / of the system, 

(4.12) 



is in general different from the bulk neutron excess 5o- Using Eqs. (|2.3| ) and (|2.4| ) the 
difference between the total volume energies reads 

\dp ddo p ) \dp d6o p J p 
1 dS 

= -?^i^~^o)A, (4.13) 

with all variables evaluated at the density Psat for ^o- The Taylor expansion of E^ c{I) about 
5o gives 

E^,,{I) = E^,M + -7^(1- ^o) + TTl^ - '^o)' + • • • • (4.14) 

p OOq 2p OOq 

Comparison of Eq. ([4.14|) with Eq. ( [4.13| ) shows that i?v,/i(5o) — E^,e{I) is of second order in 
the small quantity I — 60 and, consequently, one can approach E^^^{6q) by £'v,e(-^)- Since in 
the limit I —>■ [or 60 —>■ 0) it is E^ dl) = civ + JI^, for small 60 and / we can write 

^v,m(5o) = av + J/'. (4.15) 

Equation ( |4.13| ) is actually a surface term (i.e., proportional to A^^^). This can be seen 
as follows. Defining the effective neutron and proton volumes as [III 



N Z 

K = — and K = — , (4.16) 

Pn Pp 
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one finds 
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(4.17) 



V V A/p l-u. 

In the limit of a very large nucleus AV/V can be written as QS/V |T^, where S is the 
surface area and 9 the neutron skin thickness. With V = 4:TtR^/3 and S = AttR'^, 



AV^SQ^ 36 
V R ro 

Combining ( [4.17|) with ([4.18|) and neglecting 5^ in front of unity leads to 

30 



(4.18) 



1-5. 
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(4.19) 



which confirms the statement that ( [4.13| ) is proportional to A^^^. 

For small 5o, and recalling Eq. (|4.3| ) for G, from Eq. ( [4.19|) it is easy to show that the 
bulk neutron excess 5o and the overall neutron excess / are related by 

\ -1 / 



provided that {9J/4:Q)A'^/^ < 1. This is an interesting result, since it represents a new 



(4.20) 



manner to compute the surface-stiffness coefficient Q, in addition to Eqs. ( [4.3| ) and ( [4.6| ). 
By calculating a series of large and uncharged finite nuclei with / constant, and evaluating 
5o from the central variational densities, from Eq. ([4.20|) one can extract the value of Q for 
the considered interaction. 

Using the LDM expansion ( [4. 5] ) for i?s^(5o) plus Eq. ( [4. 20] ), and neglecting asymmetry 
effects in the curvature energy (PA^^^ terms), one recovers from Eq. ([4.11|) the liquid droplet 
mass formula for an uncharged nucleus of small asymmetry |]l|: 



E = a.A + E,,oA'/'+(^E^,o-^'jA'/' 



+ 



J 



'9J2 



2EsoL 



^-1/3 



PA. 



(4.21) 



,4g K 

In the above equation Eg^ and i?c,o are the surface and curvature energy coefficients for 
the symmetric case, and the corrective term —(2E^q/K)A^^^ accounts for the surface com- 
pression effects that appear in finite nuclei |jl]j5^. The formula ( [4.2 1|) can also be obtained 
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from Eq. ( ^.101 ) by expanding £'v,c(5o) as + J6q and using Eq. ( [4.20[ ) together with the 
expansion ( [4 .41) for i?s,e(^o)- 

Tables 5 and 6 collect the energy per nucleon, obtained in various ways, for uncharged 
large nuclei using the SkM* and NL2 interactions. The mass number ranges from A = 250 
to A = 20000, while the overall neutron excess is fixed at / = 0.2. In both tables we list 
in the second column Sq{A), which we obtain from the interior densities produced by the 
self-consistent calculation of the finite nucleus with A particles: 

= ^-W-Z^pW . (4.22) 

Pn(0)+pp(0) ^ ' 

Note from the tables that Sq approaches I a.s A grows. The energies (per nucleon) of the 
finite nuclei are given in the third column of Tables 5 and 6, labeled by EpT<{{I). In the 
next two columns we show the predictions of the mass formulae ( 4.10| ) and ( |4.11j ), labeled 



by i?MF,e(^o) and Emf,^i{So) respectively. In both cases we have included an additional term 
— {2E^q/K) to account for the compressional effect. The rightmost column £'ldm(-^) exhibits 
the output of the LDM mass formula ( [4.21|) . 

The EMF,^i{So) results agree almost perfectly with the self-consistent values Ep^{I) for 
the largest analyzed nuclei. The differences to Ep^{I) are in all cases smaller than 0.5%, a 
similar quality to that found in Ref. for a calculation performed in the symmetric case. 
Though the quahty of EMF,e{So) is also remarkably good, the overall agreement with the 
self-consistent calculation is a little worse. It was discussed in Section 3 that the surface 
tension which is minimized in the semi- infinite calculation is the Gibbs surface tension cr^, 
while (Je is not. The liquid droplet model results also agree very well with Epjsf{I), and by 
extension with -EMF,/i(^o) and Emf,c{^o)j as expected from the small value of /. This fact 
shows that Eqs. ( ^.10] ) and ( [4. Ill ) have the right limit for small 6o. 



With respect to the LDM mass formula ([4.21 ), Eqs. ( 4.10 ) and ( [4.11 ) show two main 



differences. On the one hand, the LDM volume and surface energy coefficients have been 
expanded up to quadratic terms in asymmetry, while our calculation includes it to all orders. 
On the other hand, the LDM expression neglects the dependence on asymmetry of the 
curvature energy term, which we have taken into account in Eqs. ( [4.10| ) and ( [4.11[ ). Looking 
at each coefficient separately, we have checked that the relative difference (between our 
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calculation and LDM) is much more important in the curvature energy coefficient than in 
the other ones. However, due to the A dependence of the mass formula, for very large nuclei 
most of the final discrepancy between our total energies and the LDM ones comes from the 
volume and surface terms. This is no longer the case for small nuclei, where the curvature 
term is mainly responsible for the disagreement. 

An advantageous feature of the mass formulae ( [4.10|) and (|4.11| ) with respect to the 
LDM one is that they also could be employed when drip particles exist. In such a case, one 
should use the expressions of the e- and yU-volume energies (-Ev,e and -Ev.^t) with drip. The 
compressional energy should be calculated following techniques similar to those of Ref . |61[| , 
where this correction was derived at finite temperature for helium clusters. 

The neutron skin thickness B can be extracted from the finite nuclei results through Eq. 



( |4.19| ). Reading 60 from Tables 5 and 6, we have found a good agreement with reference to 
the values of G calculated in the semi-infinite medium, specially for large A. The agreement 
worsens slightly if one defines the neutron skin of a finite nucleus to be (see Ref. for a 
discussion) 

e = i?n - i?p , (4.23) 

replacing the actual neutron and proton distributions by spheres of radii and -Rp with 
constant densities Pn(0) and Pp(0): 

^nRlp^iO) = N 

^7ri?Jpp(0) = Z. (4.24) 

It can be shown that Eq. (|4.23|) transforms into Eq. ( [4.19|) in the the limit of small asym- 
metries and to ffist order in A^^^. 

As told before, one can calculate Q from the difference I — 6q hj means of Eq. ( [4.201 ). 
From the values of Tables 5 and 6 we obtain Q = 39.5 MeV for SkM* and Q = 41.9 MeV for 
NL2. These results agree nicely with those reported in Table 2 that were obtained directly 
from the semi-infinite nuclear matter calculations, Eqs. (|4.3|) and 
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5 Summary 



In this paper we have investigated the surface properties of asymmetric semi-infinite 
nuclear matter with arbitrary neutron excess. This has been done within a semiclassical 
context, by means of density functional techniques, for non-relativistic and relativistic mod- 
els. Specifically, we have used the extended Thomas-Fermi approach including gradient 
corrections of order fi^, together with Skyrme forces in the non-relativistic case and the 
non- linear a — u model in the relativistic case. 

First, we have discussed the coexistence between a nucleus and drip particles under 
the bulk equilibrium approximation. Next, we have studied the surface properties of two- 
component systems. We have found the neutron and proton density profiles in the semi- 
infinite geometry by solving self-consistently the variational Euler-Lagrange equations, at 
given values of the bulk neutron excess. General trends of the evolution of the nuclear 
surface with the asymmetry have been obtained by exploring the surface thickness t, the 
surface width b and the neutron skin thickness O. 

We have treated the calculation of the surface and curvature energy coefficients accord- 
ing to the two definitions of the reference energy -Ercf- The self-consistent calculation of 
the density profiles corresponds to the minimization of the surface tension with the Gibbs 
prescription for -Eref- An important question is the separation of the curvature energy into 
geometrical and dynamical parts. Though in the Gibbs prescription one can avoid the dy- 
namical term by partial integrations of the Laplacian, there exists always a non-vanishing 
dynamical contribution at 6q ^ if the e-definition of i^ref is chosen. We also have calculated 
the surface stiffness coefficient Q. We have extracted it from the neutron skin B and from 
the surface energy coefficients, finding a good agreement between both methods. 

To ascertain the origin of some differences in the surface properties between the non- 
relativistic and relativistic models, we have built up a non-linear a — uj parametrization with 
the same nuclear matter properties and surface energy at zero asymmetry as SkM*. The 
non-equivalent behaviour in the evolution of the surface with asymmetry is mainly due to 
the fact that the surface-stiffness coefficient is smaller for the relativistic parametrization. 



Comparing our results for Skyrme forces with those of Kolehmainen et al. |]TT| we find 
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a good correspondence, except for the curvature energy per unit length 7e for which the 



imphcit curvature dependence was not considered in Ref. |TT[. In the relativistic model we 
could make some comparisons with the Thomas-Fermi results of Von-Eiff et al. |^|. For 
interactions like NLl and NL-SH with a small effective mass, including the Ti^ inhomogeneity 
corrections reduces the surface energy and thickness, while the surface stiffness Q increases. 

Our work has also been concerned with the analysis of the impact on the surface- 
symmetry properties of the various quantities that characterize the relativistic interaction. 
The strongest dependences have been found with the scalar mass mg, that determines the 
range of the attractive potential, and with the effective mass m^/m, that somehow reflects 
the non-local effects. 

The surface and curvature coefficients derived in this paper allow to write mass formulae 
that can be extended to the case when drip particles exist. They could be very useful in 
physical situations involving large asymmetries. We have checked the predictions of our 
mass formulae with the self-consistent energies of a calculation of uncharged large finite 
nuclei, and with the LDM mass formula in the low asymmetry limit. The agreement with 
the self-consistent calculations is as good as in the symmetric case. 

Hot asymmetric nuclear matter and finite nuclei have been investigated in non-relativistic 
and relativistic calculations. However, to our knowledge, hot semi-infinite nuclear matter has 
been analyzed only for the symmetric case in the relativistic mean field theory . A natural 
extension of our work would be to study the surface and curvature properties of asymmetric 
semi-infinite nuclear matter at finite temperature, since the combined effect of asymmetry 
and temperature has an important bearing on astrophysical objects and energetic heavy ion 
collisions. 
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6 Appendix A 



In the non-relativistic case with Skyrme forces the energy density for an uncharged 
nucleus can be written as 



1 

S = — (/nrn + /prp) + -to 



1 

16 
1 

16 



t2 1 + 



X2 



3ti 1 



1 



2 

Xi 

2 
1 



+ 12^^^" 



- -Wo (pV- J + PnV- Jn + PpV- Jp) 



(6.1) 



The quantity fq {q = n,p) is related with the effective mass of each kind of nucleon through 



_ m 
m* 



2m 88 



+ 



m 

IH n 

2n^ 



Xi 



tl 1 + -f + ^2 1 + ^ 



X2 



ti{xi + 



(6.2) 



where p = pn + Pp is the total particle density. 



In the ETF approach to order h , the kinetic energy density Lr„ and the spin density 



Jg are given by 

r 3 2 /3 /3 1 (Vp^ ^Vf Vn ^niZAl!, 

/.r, = -(Stt ) /,p/ + ^/.^^ - 3 V/,- Vp, - Y^P.-^ + 



2P.^ 



2p,, 



with 



TTl 

h = (Vp + Vp,) 



(6.3) 
(6.4) 

(6.5) 



In the relativistic formulation, the mean field Hartree energy density for an uncharged 



nucleus within the non-linear a—cu model reads [p!2| ,pmpi 



(6.6) 
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where is the third component of the isospin operator and the subindex a runs over 
occupied shell-model orbitals of the positive energy spectrum. The relativistic effective 
mass (or Dirac mass) is defined by m* = m — g^cj). Sf represents the additional contribution 
to the energy density coming from the fields 0, V and R associated with the a, u and p 
mesons respectively: 



;V0)2 + m202 



(VR? + mlR^ 



(6.7) 



In the relativistic expressions we take units h = c = 1. 

The corresponding semiclassical energy density has a similar structure to Eq. 
except that the nucleon variables are the neutron and proton densities instead of the wave 



functions. In the RETF approach it reads p8| , pT 

1 



E 



k^qe^q + kl^e^q - m*" In 



4 1 k^q + £-pq 



+ E [Biq{kFq,m*){V pqf + B2q{kFq,m*) (Vpg-Vm* 



+B-igikFg, m*){Vm*y\ + g^Vp + -gpR{p^ - Pn) + Ef 



(6.8) 



Here kpq = (SttV^)^/^ is the Fermi momentum and epg = Jk^g + m*^. The functions Big 
are the coefficients of the relativistic corrections of order h'^ to the TP model: 

r2 / 



Big{kFg,m*) 

B2q{kFg,m*) 

B3g{kFg,m*) 



TT 



24H .2 



Fq^Fq 



m 



Gk-pqS^ 



In 



epg + 2kFg In 

kpg + e-pg 



kpg + Epg 



m 



Fq 



m* 



k'^ 

'^Fq 



247r% 



- 2 + ^ 1 In 



4q\ , kpq + e-Fq 



Fq L'^Fg 



'^Fqj 



m 



(6.9) 
(6.10) 
(6.11) 



We remark that the semiclassical functionals ( |6.3|) and ( |6.8| ) do not contain any laplacian 
operators A because they have been removed by suitable partial integrations. 
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7 Appendix B 



The expression (|3.27| ) for E^l^, the dynamical part of the e-curvature energy coefficient, 
requires the evaluation of the derivatives dpg/dn (with k, the curvature). We summarize 
in this appendix how this is achieved. For further details (but restricted to the symmetric 
problem) we refer the reader to the works p[ and ||^. The starting point is the Wigner- 
Kirkwood (WK) expansion of the particle density to second order in h. In the non-relativistic 
case it reads 0,|3I| 



WK 
rg, nr 



1 /2m\2/M 

3^ yi^j \ J 



3/2 



16m 



4 n 



1/2 



-1/2 



4 n I /, . 



-3/2' 



(7.1) 



where Vq = 6T-C/6pq is the one-body potential for a Hamiltonian Ti. 
In the relativistic model we have p8 



WK 

r g, rel 



3^ 



247r2 



2£^ + 41n^^^±^l AVq 



+ 2 



k^qUl* 



^Fq 



yVa-Vm*) + 



*\2 



^Fq 



2-i|i (Vm*) 

V '^Fq/ 



2— Am* 

kFq 



(7.2) 



with k^q = (/ig — Vq)"^ — m*^ and Vq = g^V + gpT^R/2 (see Appendix A for notation). 

Since in the semi-infinite geometry ^ —>■ d/dz and A — * d'^/dz'^ + Kd/dz, one easily finds 

-1/21 



dp, 



WK 
5,nr 



Ok 



^Pq, rel 

dn 



/2my/2 



247r2 V /iV 
1 



5^ (^'^-^^ 

fq \ fq 



1/2 



/g 



127r2 



-^g , oi ^Fq + £Fq\ , ru 
-r^ + 2m — - — — - | Vq + - — m 



^Fq 



k 



Fq 



m' 



Fq 



where the primes denote a derivative with respect to z. 



(7.3) 
(7.4) 
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Following standard techniques to pass from the WK expressions to the ETF or RETF 
functional (see e.g. Refs. [p|,^^, ^|34| ) , it is possible to eliminate algebraically the deriva- 
tives of the potential Vg in favour of the derivatives of the density pg and the effective 
mass. This way, from Eqs. ( [7. 3D and ( |7.4| ) one gets the corresponding ETF and RETF 
contributions: 



dn 



ETF 



RETF 



1 

12^ 



(7.5) 



kpq + ^Fq 



k 



Fq 



+ 2 In 

\£Fq 



kpq + ^Fq 



m 



(7.6) 



Finally, replacing these results into Eq. ( |3.27| ) one is able to calculate E^^J^ for Skyrme forces 
and the a — u model. 
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Table captions 



Table 1. Properties of infinite nuclear matter for the forces used in this work. The last 
two rows show the values 5nd and 6^^ of the relative neutron excess at the neutron and 
proton drip points. 

Table 2. Properties of symmetric semi-infinite nuclear matter (surface and curvature en- 
ergy coefficients Es and E^, and surface thickness t), plus the surface-stiffness coef- 
ficient Q. The labels (1) and (2) for Q stand for Eqs. (|47^ ) and ( ^^ . Besides, we 
present results from previous works: Ref. [|l^] (KPLT) and Ref. |^ (VPSW). The 
calculations of VPSW were performed in the Thomas-Fermi approximation. 

Table 3. Change of the surface properties of symmetric matter, and of S^d and Q, with 
the properties of the relativistic interaction. The first row gives the results for the 
non-linear set NLM. This set is defined by = — 16 MeV, pnm = 0.16 fm~'^, K = 200 
MeV, m*^/m = 0.70, J = 30 MeV and = 500 MeV. (The nucleon and u- and 
p- meson masses are m = 939 MeV, rriy = 783 MeV and rrip = 763 MeV.) The next 
rows correspond to sets of parameters that differ from NLM only by the property 
listed in the first column. 

Table 4. Same as Table 3 for the surface properties of asymmetric matter at Sq = 0.212 
(below neutron drip) and 6q = 0.6 (above neutron drip). Units are MeV for the 
energies and fm for t and B. 

Table 5. Energy per nucleon of uncharged large nuclei with mass number A and overall 
neutron excess I = 0.2, for the Skyrme force SkM*. From a semiclassical calculation 
of finite nuclei we obtain the asymmetry 6q at the center of the nucleus, Eq. ( [4.22|) , 



and the energy Ep^{I). Emf,c{Sq) and -E'MF,^t(<^o) are the results of the mass formulae 
( |4.10| ) and ( [4.11[ ), including a term —{2E^q/K). Finally, i?LDM(-^) is the prediction of 



the LDM mass formula (|4.21 



Table 6. Same as Table 5 for the relativistic set NL2. 
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Table 1 





SkM* 


SIII 


NLl 


NL2 


NL-SH 


av (MeV) 


-15.77 


-15.85 


-16.42 


-17.02 


-16.35 


Pnm (fm-3) 


0.160 


0.145 


0.152 


0.146 


0.146 


K (MeV) 


216.6 


355.4 


211.1 


399.2 


355.3 


m*^/m 


0.789 


0.763 


0.573 


0.670 


0.598 


J (McV) 


30.03 


28.16 


43.46 


45.12 


36.12 


L (MeV) 


45.8 


9.9 


140.2 


133.4 


113.7 


4d 


0.322 


0.337 


0.240 


0.222 


0.284 




0.820 


0.970 


0.637 


0.749 


0.714 
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Table 2 





SkM* 


SIII 


NLl 


NL2 


NL-SH 


(MeV) 












this work 


16.00 


16.47 


17.38 


19.65 


17.22 


KPLT (1/36) 


16.05 


16.55 








KPLT (1/18) 


17.96 


18.79 








VPSW (TF) 






19.78 




20.07 


(MeV) 












this work 


10.53 


7.33 


12.59 


9.10 


8.40 


KPLT (1/36) 


10.74 


7.50 








KPLT (1/18) 


13.87 


10.10 








t (fm) 












this work 


2.23 


1.72 


2.11 


1.50 


1.50 


KPLT (1/36) 


2.26 


1.75 








KPLT (1/18) 


2.45 


1.93 








VPSW (TF) 






2.90 




2.09 


Q (MeV) 












this work (1) 


39.6 


63.8 


29.8 


41.9 


34.7 


this work (2) 


38.4 


63.3 


29.0 


41.6 


34.2 


KPLT (1/36) 


38 


61 








KPLT (1/18) 


34 


54 








VPSW (TF) 






24.4 




27.6 



43 



Table 3 



So = 







(MeV) 


Ec (MeV) 


t (fm) 


Q (MeV) 


Set NLM 


0.385 


14.81 


9.53 


1.68 


28.8 


av = -17 MeV 


0.439 


14.67 


9.26 


1.59 


30.4 


Pnm = 0.145 fm-3 


0.387 


14.29 


8.88 


1.68 


29.4 


K = 300 MeV 


0.355 


17.56 


9.61 


1.61 


30.0 


m*^/m = 0.55 


0.446 


16.54 


12.01 


2.03 


22.2 


J = 40 MeV 


0.252 


14.81 


9.53 


1.68 


35.0 


ms = 550 MeV 


0.385 


11.82 


6.05 


1.18 


38.4 



Table 4 









So = 


0.212 










So = 


0.6 








Es,e 


Ec,e 


Es,ij, 


Ec,ix 


i 


e 


Es,e 


Ec,e 


Es,iJ. 


Ecu 


t 


G 


Set NLM 


18.7 


11.6 


VIA 


8.1 


2.06 


0.43 


17.0 


18.8 


1.65 


1.00 


4.63 


0.70 


Ov = -17 MeV 


18.4 


11.2 


VIA 


8.0 


1.94 


0.40 


17.0 


18.5 


1.64 


0.95 


4.45 


0.70 


pnm = 0.145 fm" 


-3 18.1 


10.9 


11.9 


7.5 


2.06 


0.43 


16.0 


17.9 


1.58 


0.94 


4.63 


0.68 


K = 300 MeV 


21.2 


11.6 


15.0 


8.1 


1.94 


0.41 


22.6 


24.7 


2.41 


1.01 


4.17 


0.82 


m^/m = 0.55 


21.7 


15.1 


13.6 


10.6 


2.58 


0.56 


19.3 


21.4 


1.58 


1.23 


6.20 


0.91 


J = 40 MeV 


21.0 


14.4 


11.2 


6.8 


2.26 


0.51 


15.1 


19.2 


1.40 


0.99 


4.89 


0.48 


ms = 550 MeV 


14.9 


8.4 


10.0 


5.4 


1.51 


0.33 


13.8 


15.2 


1.36 


0.73 


3.69 


0.56 
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Table 5 



SkM* (/ = 0.2) 



A 


So 




^MF,c('^o) 


EMF,fi{So) 


-E'ldm(-^) 


250 


0.152 


-12.129 


-12.124 


-12.171 


-12.111 


500 


0.160 


-12.664 


-12.667 


-12.652 


-12.561 


1000 


0.168 


-13.080 


-13.085 


-13.074 


-13.068 


5000 


0.180 


-13.723 


-13.728 


-13.722 


-13.711 


10000 


0.184 


-13.905 


-13.909 


-13.904 


-13.892 


20000 


0.187 


-14.047 


-14.049 


-14.047 


-14.035 



Table 6 







NL2 


(/ = 0.2) 






A 


So 


^fn(/) 


EMF,e{So) 


EMF,fi{So) 


-E'ldm(-^) 


250 


0.137 


-12.447 


-12.456 


-12.422 


-12.530 


500 


0.147 


-13.064 


-13.079 


-13.052 


-13.111 


1000 


0.156 


-13.540 


-13.554 


-13.534 


-13.564 


5000 


0.172 


-14.276 


-14.285 


-14.275 


-14.267 


10000 


0.177 


-14.483 


-14.493 


-14.482 


-14.466 


20000 


0.181 


-14.646 


-14.659 


-14.646 


-14.624 
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Figure captions 

Figure 1. Neutron and proton chemical potentials of bulk matter (//n and /ip) as a function 
of the relative neutron excess S, for the relativistic interaction NLl. The vertical 
slashes indicate the neutron and proton drip points (^nd and (5pd) and the critical 
point {Sc) where the densities of the two phases become equal. For (5 < (5c it is 5 = Sq, 
the neutron excess of the nuclear phase. When proton drip occurs {S^, < 1) the dashed 
lines allow one to read the neutron excess of the drip phase that is in equilibrium 
with the nuclear phase at 6o. 

Figure 2. Same as Fig. 1 for the densities of the nuclear and drip phases, po and pd, against 
5^. The superimposed dot-dashed hues illustrate an example of coexistence between 
a nuclear medium with 6q = 0.6 and density po = 0.082 fm~^ and drip matter with 
Si = 0.9 and density pd = 0.064 fm"^. 

Figure 3. Equation of state for the NLl parametrization. The solid lines represent the 
pressure P as a function of the density p for several values of S (from the neutron 
drip point to neutron matter). The dot-dashed hne defines the coexistence curve. The 
long-dashed and short-dashed lines are the diffusive and mechanical stability curves. 
Some samples of points on the coexistence curve that can be in phase equilibrium have 
been joined by dashed lines. 

Figure 4. Same as Fig. 3 for the neutron chemical potential pn- 

Figure 5. Same as Fig. 3 for the proton chemical potential Pp. 

Figure 6. Schematic representation of the density profile p{z) of the semi-infinite system 
with drip particles. The quantities po and pd are the asymptotic densities when z — > 
— oo and z — > oo, respectively. 

Figure 7. Neutron and proton local density profiles of semi-infinite matter for the relativis- 
tic set NLl and the Skyrme force SkM*. They are drawn for several bulk asymmetries: 
5o = (symmetric system), Sq = 0.24 (below neutron drip). So = 0.4 (above neutron 
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drip), and finally at a 60 above proton drip. For the asymmetric cases the lower curves 
are the proton densities. The vertical bars show the surface thickness t, the distance 
where the density of the nucleus drops from 90% to 10% of its central value. 

Figure 8. Surface thickness t of the neutron (n) and proton (p) density distributions against 
the bulk neutron excess 60 for the relativistic sets NLl and NL2 and for the Skyrme 
forces SkM* and SIII. 

Figure 9. Same as Fig. 8 for the surface width b of the neutron and proton density distri- 
butions. 

Figure 10. Number of neutrons per unit area in the surface region (see text for explanation) 
as a function of 6q for the NL2 parametrization. The vertical line indicates the neutron 
drip point. 

Figure 11. Neutron skin thickness versus So- 

Figure 12. Dependence of the surface tension upon 6q for the two definitions of the refer- 
ence energy discussed in the text. 

Figure 13. Dependence of the curvature energies per unit length 7e and 7^ upon 6q for 
SkM* and SIII. The dashed lines show the dynamical contribution 7g-*''^. 

Figure 14. Same as Fig. 13 for the relativistic sets NLl and NL2. 

Figure 15. Surface tension and surface thickness for SkM* and for the relativistic set RSkl* 
adjusted to the nuclear matter properties and surface energy at 5o = of SkM*. 

Figure 16. Surface-stiffness coefficient Q against the bulk symmetry energy J for some 
relativistic and non-relativistic parametrizations. 

Figure 17. Mass-formula symmetry coefficient Osym calculated for ^°^Pb, Eq. ( [4.9|) , against 
the bulk symmetry energy J for the forces of Fig. 16. The dashed lines roughly indicate 
the empirical region. 
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